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Direct Scattering for the Benjamin-Ono Equation with 
Rational Initial Data 


By Peter D. Miller and Alfredo N. Wetzel 


We compute the scattering data of the Benjamin-Ono equation for arbi¬ 
trary rational initial conditions with simple poles. Specifically, we obtain 
explicit formulas for the dost solutions and eigenfunctions of the associ¬ 
ated spectral problem, yielding an Evans function for the eigenvalues and 
formulas for the phase constants and reflection coefficient. 


1. Introduction 


The Cauchy problem of the Benjamin-Ono (BO) equation 


du „ du 
ot ox 


dx"^ 


= 0, u(x, 0) = rio(a;), xG 


( 1 ) 


describes the weakly nonlinear evolution of one-dimensional internal grav¬ 
ity waves in a stratified fluid The variable u corresponds to the 

wave profile, e is a positive real number that quantifies the effects of dis¬ 
persion, and the operator Ti denotes the Hilbert transform defined by the 
Cauchy principal value integral 

n[u]{x,t):=-r "^df. ( 2 ) 

The Cauchy problem 0 is solvable, for suitable uq decaying as x —>■ oo, 
by an inverse-scattering transform (1ST) associated to a spectral problem 
(cf. ([^) in which uq appears as a potential. The aim of this paper is to 
construct explicitly the dost solutions and “bound state” eigenfunctions 
for this spectral problem with generic (simple poles) rational uq, and 
to deduce therefrom the corresponding scattering data. This work was 
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inspired by a brief note in a paper of Kodama, Ablowitz, and Satsuma 
[1] on the intermediate long wave (ILW) equation from which BO arises 
as a limiting case. In [1], the spectral problem © is analyzed for a 
“Lorentzian” potential (uo(x) = 2i^(l + x^)~^ for ly € R) and it is shown 
that the eigenvalues are, under some conditions on zz, the roots of certain 
Laguerre polynomials. While some complications and essential difficulties 
in extending this remarkable calculation were highlighted by Xu [5], we 
can now report that a vast generalization is indeed possible. 

It is plausible that a method capable of treating a general class of 
rational potentials can subsequently be extended to arbitrary potentials 
by an appropriate density argument. Moreover, our results allow for the 
analysis of the scattering data in the zero-dispersion limit (e —)• 0) [6]. 


2. Background 


The BO equation is the compatibility condition for a Lax pair mil], and 
a corresponding 1ST for Q was first found by Fokas and Ablowitz [9] 
with key simplifications accounting for reality of the potential given later 
by Kaup and Matsuno [10] . The first step of the 1ST solution of Q is the 
association of uq with a set of scattering data via the spectral problem 

d'w~^ 

it— - V — w~) =—uqw'^, X G M, (3) 

dx 

where A G C is a spectral parameter and the superscripts ± denote the 
boundary values of the function w{x), analytic and bounded on C\M, 
from the upper- and lower-half x-plane on the real line. 

Recall the Cauchy operators defined by the singular integrals 






'-oo y-ix± di) 


dy. 


( 4 ) 


The Cauchy operators ±C^ are bounded operators on L^(M) and are 
complementary (due to the Plemelj formula — C~ =1) orthogonal 
projections from L^(M) onto its Hardy subspaces BI^(M). For details on 
singular integrals see UDEj. Using the Plemelj formula and Liouville’s 
Theorem, ([^ is projected into the equation 

(1111 ' 

-b A - wq) = -C'^ [uow~^] , (5) 


where wq is a suitable constant |5]. Thus w~ has been eliminated, and 
([^ is related to the eigenvalue problem Cw~^ = Xw^ for the self-adjoint 
operator C := —itdx—C^UQC^ acting on BI"''(M) having essential spectrum 
and possibly some negative real point spectrum m- 
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To calculate the scattering data of one first obtains the Jost solu¬ 
tions. These are certain solutions of (also analytic and bounded for 
Im{x} > 0) that are well-defined for A > 0; in the literature [9l [TOl [TH [T5] 
they are denoted by w~^ = M(x;A) and iV(x;A) (for tco = 1) and 
w~^ = M(x;A) and N{x]X) (for wq = 0). They may be characterized 
for fixed A > 0 via their asymptotic behavior for large real x as follows: 

M —f 1, f 1 as x —)■ —oo 

- -w ( 6 ) 

The Jost solutions solve (second kind) Fredholm integral equations [9] . It 
is then possible to show that M, N, and N satisfy the scattering relation 

M{x-,X)-N{x-X) = P{X)N{x-X), AgM+, (7) 

determining a function /? : M'*' —)• C (independent of x), called the reflec¬ 
tion coefficient. 

We next introduce the “bound state” eigenfunction = ‘hj(x) G 
]HI+(M) satisfying with Wo = 0 for a given eigenvalue A = Aj < 0 and 
normalized by the condition 


x^j{x) —)• 1 as |x| —)• oo (uniformly for Im{x} > 0), 
or equivalently, as can be shown asymptotically from ([^, 


1 

2TTi 


f 

J —C 


uo(x)‘hj(x) dx 


Xj. 


( 8 ) 

(9) 


Importantly, while the Jost solutions are initially dehned only for A > 
0, M and N have analytic extensions into the complex A-plane while in 
general M and N do not [9]. Indeed, for each fixed x with Im{x} > 0, M 
and N can be shown to be the boundary values on M"*" from the upper- 
and lower-half A-planes respectively of a single function W meromorphic 
on C \ M"*". The function W has as its only singularities in A on C \ M"*“ a 
discrete set of simple poles located precisely at the eigenvalues X = Xj <0; 
see IISI Theorem 2.1]. We will henceforth use the notation IT+ (resp., 
W-) to denote the Jost solution M (resp., N). (For the inverse theory, it 
is then important that N{x\ A) may be fully eliminated from Q to yield 
a nonlocal jump condition relating the boundary values W± of W across 
its branch cut M+.) The meromorphic function W then has a Laurent 
expansion about each eigenvalue Xj of the form: 

W{x]X) = + {x + -fj) ^j{x) + 0{X - Xj) (10) 

where G C is the phase constant associated with the eigenvalue Xj. 
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Definition 1 (BO scattering data of the potential uq). The scattering 
data associated with uq are as follows. 


Reflection coefficient /3(A) for A G M"''.' defined by 0 ; see also ( pd] ). 
Eigenvalues {Aj < 0}^^; determined from the spectral problem^^) . 
Phase constants defined in terms of the eigenvalues and the 

corresponding normalized eigenfunctions by (10). 


3. Direct scattering for rational potentials 


3.1. Rational Potentials 

A real bounded potential uo that is rational with simple poles and that 
decays as x —)• ±00 necessarily has the form: 

p 

uo{x) = '^—^ -hc.c., (11) 

where {cp}^^^ are nonzero complex numbers and the poles {zp}p^i with 
Imjzp} > 0 have distinct real parts increasing with p. We impose in 
addition the condition + Cp) = 0, ensuring uq G L^(M). 

Let f{x) be the particular anti-derivative of uo{x) given by 

fix) ■='^Cp\Eog{i{x - Zp))+ ^\ +C.C., f'{x) = uo{x), (12) 

p=i f 

where Log(-) is the principal branch (|Im{Log(-)}| < tt). We denote 
by G the domain of analyticity of / as illustrated in Figure Given 
R > 0 sufficiently large, we define subdomains Gp of G n {|x| > R}, p = 
1,..., P — 1, consisting of points lying between the branch cuts emanating 
from the points Zp and Zp+i- Then we define Go as the half-plane Re{x} < 
Re{zi} and Gp as the half-plane Re{x} > Re{zp}. 


Lemma 1. For any m = 0,1,... ,P, we have 


lim /(x) = 27Ti V 

X—^00£Gm 

p=l 


Cp and 


lim fix) = —27riy^ i 




(13) 


P=1 


Proof: This follows immediately from formula (12). 


The advantage of considering potentials no of the form ( [IT| ), as first 
noted in [3], is that upon applying the Plemelj formula to the right-hand 
side of ([^, the expression C [non’"'"] can be evaluated by residues because 
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Figure 1. Left; a domain G with subdomains Gp and G* in the x-plane 
for a potential uq of the form 0 with P = 3. The zig-zagged half-lines 
denote the (logarithmic) branch cuts of f{x). Right: contours lm{x) in 
G originating at oo in the subdomains Gm and terminating at x € G. 


w~^ is analytic and bounded for Im{x} > 0. Therefore, 

p 

■ ^ CpW^ (Zp] X) 

le— -h -h = Atco + > — - - -. (14) 

ax X — Zp 

p=l ^ 

As the x-dependence on the right-hand side is explicit and elementary 
this first-order equation can be solved in closed-form, at least assuming 
w~^{zp\ A) is known for all p. 


3.2. The Jost solution M = VF+ and reflection coefficient j3 
Let A ^ 0. To obtain a formula for the Jost solution lR^(x 5 A) we take 


rco = 1 in (14), apply the boundary condition lTq_(x; A) —)• 1 as x 


—oo, 


and integrate to obtain the convergent improper integral representation; 


VL+(x;A) = --e 




r I A + J]] 

J—oo \ _ _1 


Vp{\) 


^ Z- Zp 

P=1 ^ 


dz, (15) 


where h{x; A) := Ax -|- /(x) and Vp{X) := CpW^{zp; A) for p = 1,..., P. 
Note that by Lemma h{x; A) is dominated by Ax for large |x| provided 
A 7 ^: 0. Since A > 0, the integral can therefore be made absolutely conver¬ 
gent by rotation of the contour at infinity into the lower-half x-plane. A 
similar representation can be obtained for N{x] A) = W-{x; A). 


While VL+(x;A) given by (15) satisfies the differential equation (14) 
and is analytic for all x near M, extra conditions are required to ensure 
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analyticity for all x with Im{a:} > 0. Indeed, the factor 


p 

^ih{x;X)/e ^ ^iXx/e ^^{cl-c,)l{2e) 

p=l 


appearing in (15) may have singularities for Im{x} > 0. Requiring an¬ 
alyticity of Wj^{x]X) imposes constraints on the quantities {up(A)}^^;^, 
thus eliminating them entirely from (15). 


Definition 2. imix) denotes any member of the equivalence class of 
contours in G originating at oo in the subdomain Gm and ending at the 
point X £ G; see Figure [7J 


Definition 3. denotes any member of the equivalence class of 
contours on the Riemann surface of f originating at —too in Gq with value 
/(oo) = 0 and ending at —too in Gq having encircled the singularities 
zi,... ,Zm each exactly once in the positive sense; see Figure\^ 




Figure 2. Left: contours for the configuration illustrated in Fig¬ 
ure!^ Right: contours for the same configuration. 


Proposition 1. Let A > 0. For each m = 1,... ,P, define the func¬ 
tion Hm = Hm{z] X), z £ G, by the equation 

(17) 

Z-Zpj 


{i{z - Zm))-^--'^-^iim{z- A) := [ A + 

\ p=i 


(Note that Hm is analytic not only on G, but also in a neighborhood of 
the point Zm-) Also, let icm = e(/tm + ajm), where pm = [Re{tCm/e}] 
denotes the integer (floor) part ofKe{icm/e} and the remainder satisfies 
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0 < Re{a;m} < 1- Then, the function W'+(x;A) given by (15) is analytic 
in the upper-half x-plane if and only if for each m = 1,..., P, either 





0 when Re{icm/e} < 0 , (18) 



dz = 0 

when Re{icm/e} > 0 


and ojm / 0, (19) 


or 

Res - -= 0 when Re{icm/e} > 0 and ujm = 0, (20) 

Z=z^ [z - Zmr^^^ 

where io{zm) (see Definition^ originates in at —ioo in Gq. 

Proposition is proved in the Appendix. Note that the singularity at 
z = Zm arising from the exponential on the right-hand side of 

( [ItI ) is explicitly cancelled from Hm{z', A) by the factor {i{z — 
on the left-hand side. The essence of the proof is to deform the contour 
of integration in ( |15[ ) to pass through the point Zm when the integrand 
is integrable there, i.e., Re{zcm/e} < 0. Thus, we deal separately with 
integrals on io{zm) and on the line from Zm to x. If the integrand is 
not integrable at Zm, we first integrate by parts repeatedly to achieve 
integrability, or invoke the Residue Theorem if the singularity is a pole. 

It is useful to rewrite the conditions for analyticity of lT+(x; A) given 
in Proposition by eliminating the derivatives of Hm as follows. 


Corollary 1. Let A > 0. The function lT+(x; A) defined by (15) is 
analytic in the upper-half x-plane if and only if {vp{X)}p^i satisfy 




'C> 


E 

vP=l 


Up (A) 


Z - z„ 


-b A dz = 0, m = l,...,P, 


( 21 ) 


where Off, := unless icm!^ is a strictly negative integer, in which 
case := io{zm)- (The contours extend to —ioo in Gq making the 
integrals absolutely convergent.) 


Proof: Let Lm be a contour on the Riemann surface of / beginning 
(with value /(oo) = 0) and ending in at —ioo in Gq, and encircling only 
the singularity Zm once in the positive sense. Now, in the exceptional case 
when = io{zm), © is just the condition ( |18[ ) rewritten using ©• 
In (18) with Um 7 ^ 0 as well as in (19), the integrand has an integrable 
algebraic branch point at z = Zm so up to a nonzero factor of the form 
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(1 — a)/2 the integrals appearing in (18) and (19) can both be rewritten 
equivalently with the contour £o{^m) replaced by Lm- In the case of 
(19) we then integrate by parts + 1 times on Lm, which produces no 


boundary terms for A > 0 due to exponential decay of the integrand at 
—ioo. Eliminating Hm by means of (0 then gives ( [^ with replaced 
by Lm in these cases. The final condition of Proposition ( |20[ ), can of 
course be written as an integral about a small loop centered at z = Zm, 
and applying Cauchy’s Theorem under the condition A > 0 allows this 
loop to be deformed into the contour Lm (in the case of a locally single¬ 


valued integrand), so (20) also takes the form (21) with Cm replaced by 


Lm after elimination of Hm{z',\). Finally, the replacement of Lm with 


Um in each case to arrive at (21) amounts to multiplication of the system 


(21) by a triangular nonsingular matrix. 


Note that the conditions for analyticity of W+{x\X) as expressed in 
Corollary take the form of a square linear system of equations on the 


unknowns {up(A)}p^]^ assembled in a vector v(A) := (ui(A),... ,vp{\)) 


T. 


A>(A)v(A) = b>(A), A>0, 

where A^(A) G and b^(A) G have components 

^-ih{z;\)le 


( 22 ) 


A> 

■^mp 


(A) := / 

Jc 


dz and h^{\) \=—X / e dz. (23) 

c> z, — Zp Ju> 


The reflection coefficient /3(A) is now easily calculated. Combining the 
scattering relation ([^ with the boundary conditions ([^ gives 

ld{X) = lim {W+{x- A) - 1) for A > 0. (24) 

X^OO 

For a rational potential uq of the form 0, we may explicitly write 


/3(A) = -e 




f 


p=l-p / g 

—oo 


—ih{z\\)/e 


Mz) - 


p=i 


Up (A) 


Z- Zr, 


dz, (25) 


using equation (15) for the Jost solution VF+(x;A) and Lemma [Tj The 

.-. L-l 


coefficients {vp{X)}p^i are determined by the linear system 


3.3. Eigenfunctions and eigenvalues 

Eigenfunctions w~^ = <h(x; A) corresponding to eigenvalues A < 0 are 
nontrivial solutions in BI^(M) of ([^, or equivalently (14) for uq of the 
form 0, with ruo = 0. Integrating (|14[) with tco = 0 and imposing 
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A) —f 0 as X —)■ —oo shows that <^(x; A) necessarily has the form 


A) = —e 


_ ih{x\\)le 


^-ih(z\X) I e 


E 

p=i 




Z- Zr, 


dz, 


(26) 


where 4>p{X) '■= Cp<h(zp;A) for p = Since A < 0, the conver¬ 


gent improper integral on the right-hand side of (26) becomes absolutely 


convergent if the contour is rotated at infinity into the upper-half x-plane. 


While (26) satisfies (14) and is analytic for x in a neighborhood of 


M, to obtain the required analyticity for all x with Im{x} > 0 requires 
conditions on both A < 0 and {(?!)p(A)}^^Q. As we will see these conditions 
will also guarantee that 4>(-; A) G i.e., that it vanishes as |x| 

in all directions of the upper half-plane. 


oo 


Definition 4. denotes any member of the equivalence class of 
contours in G originating at ioo in Gq and ending at ioo in Gmi see 
Figures and 

Proposition 2. Let A < 0. The function <h(x;A) is analytic in the 
upper-half X-plane if and only if satisfy 


L 


! e 


E 

p=i 


4>p{X) 


Z - Zr, 


dz = 0, m = 1,..., P, 


(27) 


where G^ := unless icm!^ is a strictly negative integer, in which case 
Gm •= f’o(^m)- (The contours Gff extend to ioo in Gq and Gm making 
the integrals absolutely convergent.) 

The proof of Proposition is similar to that of Proposition and 
Corollary and is omitted for brevity. Introducing the vector <MX) = 
{4>i{X), ..., (l)p{X))'^ of unknowns, the conditions of Proposition take 
the form of a square homogeneous linear system A^(A)0(A) = 0 where 
the matrix A^(A) G has components 


4< 

■^mp 


(A) := / 

JC< 


^-ih(z\\)/t 


dz, A < 0. 


Z- Zr, 


(28) 


Corollary 2. The function 4>(x;A) given by (26) is both nontriv¬ 
ial and analytic in the closed upper-half x-plane if and only if D{X) := 
det(A<(A)) = 0 and <^(A) is a nontrivial nullvector of A‘^{X), unique up 
to a constant multiple. 

Proof: This follows from Proposition It only remains to check that 
dim (ker (A^(A))) = 1 whenever det(A'^(A)) = 0, a fact shown in [T7j. ■ 
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Corollary 3. Suppose D{X) = 0, and that the components 

of the nullvector 4>{X) are normalized to satisfy 


'^(fp{X) = A < 0, 

p=i 


(29) 


then the function $(x;A) satisfies the normalization condition 
$(-;A) G H+, and therefore <I>(x;A) is an eigenfunction of equation (15). 

Proof: By Corollary]^ the function <h(x;A) is analytic in the closed 
upper-half x-plane. From ( |26[ ) one then invokes the Residue Theorem to 
see that ( [^ implies (§. The fact that 4>{X) can be scaled to satisfy 
i.e., that ker(A''(A)) is not orthogonal to (1,1,..., 1)''', is proven in [T7j. 
Lastly, to show that G BI+, we check that 0 as |x| —)• oo in the 


closed upper-half x-plane. First, applying Jordan’s Lemma to (26) shows 
0 when x 


that d> 


oo anywhere in Gq. To let x 


oo elsewhere in G, 

we rewrite the limiting contour integral in (26) as lima,^ooeG™ /_oo ~ Iu< 
whenever icm!^ is not a negative integer. From Proposition this limit 
is identically zero for all m = 1,..., P, since each integral on vanishes 
when D{X) = 0. ■ 


Corollary 4. Let D{X) = 0 for some X = Xj < 0. The eigenfunction 
$j(x) := <l)(x;Aj) given by (26) with the conditions from Corollary^ 
satisfies the asymptotic condition (Isj). 


Proof: Applying I’Hopital’s rule to (26) and using h'(x; A) = A-|-tto(3:), 


lim x<h(x; A) = lim 

X^OO X^OO 


. 9 

IX 

e -p iXx + ixuo{x) 


p 


E 


</-p(A) 


i^r/)p(A), (30) 


and the desired result then follows from ( |29[ ). ■ 

The stronger result <l)(x; A) = x~^ + O as |x| —>■ oo (uniformly 

in the closed upper-half x-plane) is proven in [T7j . 

Corollaries [2||^ show that the equation D{X) = 0 involving the deter¬ 
minant of the PxP matrix A'' (A) is exactly the condition that A = Aj < 0 
is an eigenvalue for the rational potential uq of the form ( |11[ ), and that 
‘hj(x) := <l>(x;Aj) is the corresponding normalized eigenfunction. 

Remark 1. If the coefficient icm/e is a positive integer N^, the inte¬ 
gral ( [^ can be calculated by residues, and thus A<p(A) = e“*^^™/'^Pp(A), 
where Pp(A) is a polynomial in A of order at most Nm- If this holds for 
all m = 1,..., P, then P(A) is proportional via a nonvanishing exponen¬ 
tial factor to a polynomial of order at most nm=i If P = Ij this 
polynomial is a scaled Laguerre polynomial; see [1]. 









Direct Scattering for the Benjamin-Ono Equation with Rational Initial Data 


11 


3.4- Analytic continuation, Evans function, and phase constants 


The matrix (A) defined for A < 0 by (28) has an analytic continuation 


to the maximal domain C \ M"*". Recalling that h(z;X) is dominated for 
large | 2 :| by the term Xz, it is clear that this continuation is afforded simply 
by rotating the infinite “tails” of the integration contours so that they 
tend to complex oo in the direction arg( 2 ;) = — arg(fA). In the process 
of rotating the contours from their initially upward vertical configuration 
when A < 0, the function f{z) appearing in the integrand via h{z;X) 
must be analytically continued through its vertical branch cuts as well. 
Therefore we observe that D{X) := det(A^(A)) is an Evans function for 
the spectral problem with rational potential uq of the form 0 ; it is an 
analytic function in the domain C\M+ (complementary to the continuous 
spectrum for the problem) whose roots are precisely the eigenvalues. 

If the analytic continuation of A^(A) is carried out through the upper- 
half A-plane to the cut M"'', the contour C))) will be rotated through the 
left-half 2 :-plane to coincide precisely with C^; compare the two panels 
of Figure]^ Thus, the analytic continuation of the matrix A^(A) from 
A < 0 through the upper half-plane to A > 0 coincides with the matrix 
A^(A) defined in (23). In a similar way, the vector b^(A) defined for 
A < 0 with components 


b<{X):=-x[ dz, A<0, (31) 


has an analytic continuation to C \ M"*", taking a boundary value on M"*" 
from the upper half-plane that coincides with the vector b^(A) also de¬ 
fined in (23). This shows that indeed for each x with Im{x} > 0, the Jost 
solution IT+(x; A) given by the formula is the boundary value on M"*" 
from Im{A} > 0 of a function W{x; A) analytic for A G C \ M"*" with the 
possible exception of the eigenvalues A < 0 satisfying D{X) = 0; only at 
these points do the quantities {up(A)}^^]^ entering into (15) become inde¬ 


terminate as the matrix of the analytic continuation of the system 
through 0 < arg(A) < 27r becomes singular. This idea allows us to de¬ 
duce the remaining scattering data corresponding to uq of the form 0, 
namely the phase constants corresponding to the negative eigen¬ 

values {Aj}j^^. The first step is to determine how analytic continuation 
of the system (22) fails near an eigenvalue X = Xj. 


Proposition 3. Let v(A) be the unique solution of A^(A)v(A) = 
b^(A) for each A G C \ M"'' for which D{X) 0. Then each compo¬ 
nent Vp{X) is analytic in C\ M"'' except at the eigenvalues which 
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are simple poles, with corresponding Laurent expansion 

~ ^ + O (A — Aj) as A —> Xj (32) 

A — Xj 

where Fj is a constant (independent of p) given by 


r 

'■“ 2A 




^ A (\ \ iem'^b<'(A 

2 _)^. Z_^ ^p^pi^j) 




^ p=i 


2 m'''b<(Aj) ’ 


(33) 


and where {(i>p{Xj)'\p^i are the components of the right nullvector of 
A“^(Aj) normalized by ( [2^ and m is a nonzero left nullvector of A"'{Xj). 

We include the proof of Proposition in the Appendix. The following 
result then shows that the phase constant associated with each eigenvalue 
Xj < 0 by means of the Laurent expansion is given by 7j = Pj. 


Corollary 5. The function W{x-,X) whose boundary value from 
Im{A} > 0 on M+ is the Jost solution lL + (x; A) is meromorphic on C\M+ 
with simple poles only at the eigenvalues at each of which it has 

a Laurent expansion of the form (10) with <hj(x) := <h(x; Aj) defined as 


in Section 3.3 and with 7 j = Pj defined as in (33). 


Proof: The function VP(x;A) is given by (15) with the coefficients 
{i;p(A)}^^;^ obtained by solving A^(A)v(A) = b^(A), and with the con¬ 
tour of integration rotated appropriately to ensure absolute convergence 
of the integral. Substitution of (32) into ( [T^ and using (26) yields, after 
some manipulation, (10) in which 7^ = Pj- ■ 
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Appendix: Proofs 


Proof of Proposition^ Fix A > 0. The conditions of the proposition are 
equivalent to VF+(x) = VF+(x;A) being analytic at each possible singu¬ 
larity zi,... ,zp in the upper-half x-plane. For any m = 1,..., P, we use 
defining Hm{x) = Lhmix', A) in ( [T^ to express Wp{x) in the form 

lF+(x) = gm{x){i{x - 


[ {i{z-Zm)) ^Hm{z)dz, (Al) 
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where we have also introduced a function gm{x) analytic and nonvanishing 
a,i X = Zm dehned precisely by the identity egm{x){i{x — = 

iQ^h{x)/t^ We will now use Cauchy’s Theorem to deform the path lo{x) of 
integration so as to pass directly through the branch point Zm on the way 
to a nearby point x G G. While Hm{z) is analytic at z = Zm, the other 
factor in the integrand of (Al) may not be integrable at z = Zm- 

First suppose that Re{icm/e} < 0. The factor {i{z — is 

integrable a.i z = Zm so we can apply Cauchy’s Theorem to obtain 


W+{x) = gm{x){i{x - [ {i{z - Zm)) ^ Hm{z) dz 

+ gm{x)ii{x - ZmW'"^’' [ ii{z - Zm))~"'''"^''~^Hm{z) dz. (A2) 

The second term on the right-hand side is analytic for x near Zm- Indeed, 
substituting for Hm{z) its power series in i{z — Zm), 

OO 

ddm{z) = ^ ^ Vmk {'^{z Zm)) j (-^ 3 ) 

k=0 

for small enough \z — Zm\-, term-by-term integration yields 


gmix){i{x - J (i(z - Zm)) ^ffm(z) dz 


-igm{x) ^ 


k=0 


k - iCm/e 


{i{x-Zm))’^- (A4) 


Now, (A4) is analytic ai x = Zm as it is a product of an analytic function 


and a convergent power series. So, requiring that W+(x) be analytic at 
X = Zm is equivalent to demanding that the first term on the right-hand 
side of (A2) be analytic at x = Zm- But this term is the product of (i) 
gm{x) nonvanishing at x = Zm, (h) {i{x — ZmSf^'^^'^ which blows up as 
X ^ Zm since Re{zcm/e} < 0, and (hi) the x-independent integral on 
Hence, the analyticity of IF+(x) at x = Zm is equivalent to the 
vanishing of this integral (cf. &■ 

Next let Re{fcm/e} > 0 with Um / 0. Since (z(z — Zm)) ^ is 

branched and nonintegrable at z = Zm, we integrate by parts fim + 1 
times (note gm > 0) to restore integrability before deforming the contour 
ioix). With Km ■= f^’"+^r(wm)/r(wm -R fim + 1) 7^ 0 we have 


(i(z - - 


Zm)) 


—cu„ 


(A5) 
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Integrating (Al) by parts using (A5), note that {i{z — Zm)) is inte- 


grable at 2; = Zm (0 < Rejwm} < 1), so Cauchy’s Theorem gives 




W+{x) = {-l)^"'"^^Kmgmix){i{x - 




' Zm Zm)) 

1 


dz 


+ {-l)^^-+^Kmgm{x){i{x - Zm))'^-!^ P 

J Zrr 

+ Kmgm{x){i[x - Zm)y'^’"^'''^i-'^)^H^\x)-^^^^_^ _ Zm))*" 


A:=0 


(A6) 


Now, recalling icm!^ = gm + zjm-, note that the terms on the third line 
of (A6) are analytic at x = Zm- Likewise, with the help of the uniformly 
convergent power series (A3) and term-by-term differentiation and inte¬ 
gration (assuming only that \x — Zm\ is sufficiently small), the second line 
of (A6) is analytic at x = Zm- The first line of the right-hand side of 
(A6) is the product of (i) Kmgm{x) analytic and nonzero at x = Zm, (ii) a 
branched factor (i(x — Zm)Y^’^^’^, and (iii) the definite integral on ^0(2:™.)- 
This term of W+(x) is therefore analytic if and only if the integral on 
^o{zm) vanishes (cf. (EH))- 

If Re{icm/e} > 0 but instead Um = 0, then (Al) becomes 


^^+ix) — 'igm{x)(x Zm)^ 


Hm{z) 


it.o{x) (z 


dz, Hm — 0,1, 2, 3,.... 

(A7) 


Since gm and Hm are analytic at Zm, IT+(x) will be analytic at Zm exactly 
when the integral factor is single-valued as a function of x (cf. (20)). ■ 


Proof of Proposition^ Since A(A) = A^(A) given by (23) and b(A) = 
b^(A) given by ([M|) are analytic for A G C \ M"'', the solution v(A) of the 


system A(A)v(A) = b(A) will be analytic in the same domain with the 
possible exception of the eigenvalues at which D{\) = 0. In [T7] it is 
shown that these isolated singular points are simple poles of v(A), which 
therefore has a Laurent expansion about A = \j of the general form 


v(A) = 


vt-^1 

j 

X-Xi 


+ v^- ^ -p (A — Aj)vj ^ -p O ((A — Aj)^) as A —> Xj. (A8) 
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Substitution into A(A)v(A) = b(A) implies: 

= 0 , 

A(Aj)vj°' = b(Aj) - A'(Aj)v{“^', 


A(Aj)vW = b'(Aj) - ^A"(Aj)vJ. - A'(Aj)vj° 


(A9) 

(AlO) 

(All) 


Since A(Aj) is singular and has maximal rank P — 1 (see |17j). the first 
two equations (A9) and (AlO) have the general solutions 

= a[“^^0(Aj) and = p + a^°'0(Aj), (A12) 


where p is a particular solution of (AlO) and 0(Aj) is the nullvector of 
A(Aj) normalized according to (29). Here o;^- and are constants 
determined by the fact that the right-hand sides of the equations (AlO) 
and (All) lie in the range of the singular matrix A(Aj). Therefore, mul¬ 
tiplying these latter two equations on the left by m"'' where m is any 
nonzero left nullvector of A(Aj) (unique up to scaling) and using (A12) 
gives 

m'''b(Aj) 


= 

^ uPA>{Xj)(t){\j] 


and 


= 


(h'{Xj) - ^^A"{Xj)(t){Xj) - A'(Aj)p^ 


A'{Xj)4>{Xj 


(A13) 


(A14) 


To simplify these formulas, we differentiate (|28|) to obtain the identity 

(A15) 


A'(A) = ^B(A) - ^A(A)Z, 


where B(A) := [b(A) • • • b(A)] and Z := diag(zi,..., zat). Since 
m'''A(Aj) = O"'', (A13) becomes 


[-1] Ajm'''b(Aj) . m'''b(A 

“i *^m''’B(Aj)0(Aj) *^m''’b(A 


jj 


- = —ie, (A16) 

'KjJ 111 uy/VjJ 

due to (29). (That m'''b(Aj) 0 is proven in [17] .) Now, substituting 

(A12), (A15), and ( A16| ) into (AlO) and using ( [^ again yields 

A(A,)p = A(A,)Z(/,(A,), (A17) 


and hence the particular solution p may be taken as p = Z(^(Aj). Finally, 


substituting (|A15) and its derivative into (A14) shows that = Tj as 


defined by 
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